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Abstract 

This  paper  treats  the  robust  detection  problem  for  the  case  where  the  noise 
mixture  model  is  not  strongly  unimodal;  that  is,  has  heavier  than  exponential 
tails.  Under  a  constraint  on  the  support  of  the  contaminating  density,  an 
optimal  detector  is  derived.  Its  performance  is  investigated  using  both 
computer-generated  noises  and  noise  records  from  underwater  noise  fields. 

1  -  Introduction 

Robust  detection  and  estimation  techniques  have  been  developed  for  environments  in 
which  there  is  some  uncertainty  as  to  the  exact  distribution  of  the  noise  processes 
present[l,2,3].  In  [4,5]  the  problem  of  detecting  a  signal  in  contaminated  noise  is  discussed 
and  the  robust  detector  in  terms  of  risk  is  found  to  be  a  censored  probability  likelihood  test. 
The  robust  location  estimator  and  small  signal  detector  for  nominal  noise  distributions  hav¬ 
ing  lighter  than  exponential  tails  (corresponding  to  strongly  unimodal  densities)  are  described 
in  [3]  and  [6],  respectively. 

Huber  in  [3]  developed  the  basic  theory  for  robust  estimation  of  location  for  convex 
classes  of  distribution  functions.  The  motivation  for  this  work  was  the  the  desire  to  find  a 
technique  which  diminished  the  inordinate  influence  of  outliers  on  estimates.  He  found  that 
this  could  be  done  by  using  the  system  which  would  give  the  MLE  for  location  if  the  noise 
had  the  distribution  with  the  smallest  Fisher’s  Information  in  the  given  convex  class  of  distri¬ 
butions.  The  discussion  in  §2  lifts  his  restriction  on  the  distributions;  it  contains  an  investiga¬ 
tion  of  the  analogous  robust  detector  and  location  estimator  for  the  case  where  the  nominal 
has  heavier  than  exponential  tails.  The  paper  concludes  with  the  results  of  some  simulations 
which  were  performed  with  the  robust  detector  developed  in  §2. 

We  will  discuss  robust  detection  of  a  small  signal  in  independent  identically  distributed 
mixture  noise.  The  distribution  /x(x)  of  the  noise  is  given  by 

Mx)  =  n^(a;,)  =  ni(l-c)G'(a^.)  +  0<€<1  (1.1) 

i-l  i-1 

where  G  is  the  known  nominal  distribution,  H  is  the  unknown  contaminating  distribution 
and  F  is  the  actual  distribution  of  the  individual  observations.  The  univariate  mixture 
model  is  given  by 


F  =  (!-€)(?  +  ci/  , 

0<£<1 

(1.2) 

where  G  is  the  known  nominal  distribution  and  H  is  the  unknown 
more,  if  G  and  II  have  densities  g  and  h ,  respectively,  we  can  write 

contaminant. 

Further- 

/  =  (l-€)ff  -f  eh  , 

where  /  is  the  density  of  F. 

0<€<1 

(1.3) 

<!>■ 

'  Cs 

Z 
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2  -  Robust  Estimation  and  Detection  for  Noise  with  a  Non-Strongly 

Unimodal  Nominal  Density 

In  [3,4]  the  least  favorable  dbtribution  was  developed  for  the  mixture  class  of  distribu¬ 
tions  which  have  a  strongly  unimodal  nominal  density.  This  type  of  density  is  limited  to 
those  which  have  lighter  than  exponential  tails.  However,  many  of  the  noise  environments 
commonly  encountered  in  atmospheric  and  oceanographic  applications  have  been  found 
[7,9,10]  to  have  distributions  which  are  much  heavier  tailed  than  exponential.  It  therefore 
seems  proper  to  investigate  the  robust  problem  for  very  heavy  tailed  distributions. 

We  will  find  the  robust  small-signal  detector  and  M-estimator  for  the  mixture  class  of 
distributions  having  what  will  be  called  non-strongly  unimodal  nominal  densities.  This  class 
b  defined  as  follows. 

Definition  1.  A  unimodal  density  g  is  called  non-strongly  unimodal  if  it  satisfies  the  follow¬ 
ing  conditions:  ^ 

i)  sjrn[ — - — (a:— c)]  =  s^n[a:— c]  some  real  c. 
g 

a)  There  exist  Pi  and  with  P^  <  c  <  Pi  such  that 

a)  For  X  <  P^  and  x  >  Pi, then  —g’{x)/g{x)  is  monotonically  de¬ 
creasing. 

b)  For  P2  <x<  Pi, then  —g'{x)/g{x)  is  monotonically  increasing. 

Hi)  lim  g'{x)/g{x)  =  0. 

I  *  !-*<» 

For  simplicity  we  will  assume  throughout  this  section  that  c=0.  Note  that  if  ^  b  symmetric 
and  c=0  then  P^  =  —Pi. 

Cauchy  nobe  is  an  example  of  a  dbtribution  having  a  non-strongly  unimodal  density. 
Unfortunately,  the  moments  of  the  Cauchy  distribution  are  undefined.  However,  during  many 
practical  applications,  the  extreme  taib  would  be  truncated  by  the  signal  processor  so  that  we 
can,  if  it  seems  appropriate,  model  a  noise  process  as  Cauchy  within  some  finite  interval  of 
the  observation  space.  Over  such  an  interval  all  moments  exist  and  are  finite. 

The  members  of  the  Johnson  family  [7,9,11]  of  dbtributions  also  have  non-strongly  uni¬ 
modal  distributions.  Willett  [9]  found  that  a  member  of  this  class  closely  modeled  the  noise 
produced  by  shrimp  along  the  ocean  floor.  Johnson  nobe  is  generated  by  the  following 
memoryless  transformation  on  Gaussian  noise: 


I  =  X  sinh 


Here  2  b  an  observation  from  a  unit  normal  process  and  x  is  an  observation  from  the  gen¬ 
erated  Johnson  process.  The  parameter  6  controb  the  tail  weight;  very  large  6  produces  a 
nearly  Gaussian  process;  and  X  is  a  scale  parameter  given  by 


X  = 


2o^ 


_2_ 


-  1 


1/2 


where  b  the  variance  of  the  Johnson  process.  Because  6,  and  hence  the  tail  mass,  b 
allowed  to  vary  continuously,  the  Johnson  family  can  afford  a  good  basb  for  approximating 
distributions  with  very  heavy  tails. 

We  will  now  investigate  the  properties  of  the  absolutely  continuous  mixture  density 
function  /  _ 

/  =  (!-€)</  +  f*  (2.1)  ~ 

where  <7  is  a  non-strongly  unimodal  density  and  h  b  any  density  function  for  which  the 

resulting  /  is  absolutely  continuous.  We  will  also  assume  in  this  section  that  /  b  finite,  uni- _ 

modal  and  attains  its  maximum  at  1=0.  This  assumption  b  reasonable  since  many  distribu-  ’.odes 
tions  encountered  in  nature  are  unimodal  and  we  place  the  mode  at  zero  without  loss  of  gen-  — • 
erality.  Hence,  we  have  added  the  constraint  that  €h'{x)  is  less^gr^ter)  than  — (1— €)^'(a:)  for 


□  Q 


X  greater(less)  than  zero. 

Fact  1.  The  density  f  is  not  strongly  unimodal. 


Proof.  From  the  definition  of  non-strongly  unimodal  we  have 


lim 

I  1— >00 


^  =  0. 


Then  for  any  6>0  there  exists  finite  M>0  such  that  —g'{x)/g{x)  <  6  for  all  x>M .  Pick 
zq>M .  Since  both  /  and  g  are  finite  we  have  —log  /  (20)  +  log  <7(^0)  =  r  <00.  Furthermore, 
by  the  positivity  of  h,  —log  f  (x)  <  —log  g{x)+c,  for  all  x,  where  c=  —log  (1— c).  Assume  / 
is  strongly  unimodal;  that  is,  —log  /  is  convex  and,  for  positive  x,  dominates  some  line  with 
constant  positive  slope  k.  Pick  any  6  between  zero  and  k.  Then  for  Zi>Zq 


-log  /(zi)  4-  log  /  (20)  >  k{zi-Zo) 


(2.2) 


and 


-log  gizi)  4-  log  g{zo)  <  5  (21-20).  •  (2.3) 

Combining  (2.2)  and  (2.3)  yields 

-log  /  (21)  4-  log  5(21)  >  -log  /  (20)  4-  log  g{zo)  +  {k-S){zi-Zo)  (2.4) 

=  r  4-  {k-S){zi-Zo). 

But  the  right  side  of  (2.4)  can  be  made  as  large  as  desired  by  choosing  21;  hence  it  can  be 
made  larger  than  c  and  the  positivity  condition  for  h  will  be  violated  which,  by  contradic¬ 
tion,  completes  the  proof. 

Fact  2.  If  —f  /f  has  a  limitf  then 

lim  ^  =  0  (2.5) 

I  *  !-►'»  / 


The  proof  of  Fact  2  is  virtually  identical  to  that  of  Fact  1. 

In  light  of  Facts  1  and  2  we  see  that  we  cannot  contaminate  a  non-strongly  unimodal 
density  and  produce  one  which  is  strongly  unimodal.  Furthermore,  if  /  is  not  strongly  uni¬ 
modal  then  the  derivative  of  — /'//  must  be  negative  for  some  positive  arguments.  In  addi¬ 
tion,  since  /  is  unimodal,  —f'{x)/f{x)  is  zero  at  the  origin  and  positive  for  all  positive 
values  of  i;  therefore  its  derivative  cannot  be  negative  for  all  positive  x.  As  a  result  it  is  not 
possible  to  fulfill  the  conditions  of  the  lemma  of  the  previous  section  if  the  nominal  density  of 
/  is  non-strongly  unimodal. 

Fact  3.  If  the  nonlinearity  ri[x)  is  zero  at  x=±oo  and  at  x=0  and  if  sgn  [r/(a:)]  =  sgn  [i]  then 
Eprf  >0. 

Proof  Since  r/(i)>0  for  i>0  and  r}{x)<0  for  i<0,  then  rf  {x)  must  be  positive  for  some 
interval  {b,a)  about  the  origin.  Furthermore, 

00 

/>/  =0 

—00 

Then,  after  separating  the  real  line  into  (— oo,6j,  (6,0)  and  (0,00)  we  have 

a  i  00 

=-/»/  -/»/•  (2-6) 

i  —00  a 

By  the  unimodality  of  /  we  have 

00  Oh  a  00 

/-)'/  >f(l>njv'  +  /-/l+/(«)l/>/  +/'/l=0  (2.7) 

—00  h  —00  0  a 

t  For  »n  example  of  a  mixture  density  /  for  which  — / /  has  no  limit  see  [12]. 
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which  completes  the  proof. 

Let  us  now  apply  Fact  3.  We  are  given  two  distribution  functions  Fj  and  both  of 
the  form 

F{  =  (l — i)G  +  iffi,  1=1,2. 

with  density  functions  fsubt{x)  —  dFi{x)/dx.  Let  A  =  {x:  rf{x)<Q)  and  suppose 

Jhi  =  1  and  Jh^  =  0. 

A  A 


Then 


Fjr^r/  <  Ep^ri' . 

where  as  a  consequence  of  F act  3 

\EFr/\=Eprf,  i=l,2. 


(2.8) 


Hence  the  numerator  of  the  expression  for  efficacy  and  the  denominator  of  that  for  asymp¬ 
totic  variance  (see  [3,12])  are  minimized  by  F  ^  Similarly,  if  we  let 
B  =  {x:  |»y(a:)I=  a  maximum}  and  suppose 

Jhi  =  l  and  Jh2  =  0. 

B  B 


Then 

<  Erjf.  (2.9) 

Therefore,  the  denominator  of  the  expression  for  efficacy  and  the  numerator  of  that  for 
asymptotic  variance  are  maximized  by  Fj. 

If  A  =  B  then  r)  must  be  decreasing  on  the  set  where  it  is  at  its  maximum,  a  condition 
which  can  only  be  satisfied  if  the  A  and  B  are  identical  sets  of  discrete  points.  As  a  result, 
the  contaminant  would  put  all  of  its  mass  on  a  set  of  Lebesgue  measure  zero  and  would 
therefore  be  degenerate.  Since  we  have  constrained  the  problem  so  that  a  degenerate  contam¬ 
inant  is  not  allowable  we  can  state  Fact  4  as  follows. 

Fact  4.  Let  C  jbe  the  set  of  all  mixture  distributions  having  finite,  absolutely  continuous,  uni- 
modal  densities.  If  r}  satisfies  the  conditions  given  in  Fact  S,  then,  for  a  non-strongly  unimodal 
nominal  density,  it  is  not  possible  to  maximize  Eptfi'  and  minimize  Eprf  with  the  same  distri¬ 
bution  F^  C . 

The  underlying  problem  that  Fact  4  formalizes  is  that  we  have  no  clear  cut  method  by 
which  to  decide  whether  to  put  the  contaminant  mass  where  the  nonlinearity  is  large  in  mag¬ 
nitude  or  where  its  derivative  is  negative.  We  will,  in  the  remainder  of  this  section,  investi¬ 
gate  the  robust  problem  under  the  constraint  that  the  support  of  the  contaminant  is  limited 
to  an  interval.  By  limiting  the  allowable  class  of  contaminants  in  this  fashion  we  will  be  able 
to  find  a  nonlinearity  ^  distribution  Fq  =  (1— e)G  -f  cHo  which  satisfy 

min  Eprjff  =EF^rjo,  (2.10) 

r  €  V# 

max  Ef  r}^  =  Ep^Vo  (2.11) 


and 


Vo  = 


where  C*  contains  all  mixtures  with  contaminants  whose  support  is  confined  within  some 
specified  finite  interval.  Clearly,  if  Eqs.  (2.10)  and  (2.11)  are  satisfied,  then 


min/(F)  =  /(Fo), 

r  t'  # 


that  is,  Fq  is  the  minimum  Fisher’s  Information  distribution  in  C0. 
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The  robust  detector  and  estimator  for  the  class  C* 

Let  be  a  non-strongly  unitnodal  density  function  as  in  Definition  1  with  c^.  Define 
the  function 


P{x)  =  1 


g{ai)e 

gix), 


ai<z<di 

02>a:  >^2 


otherwise 


where 

k  =  -=^(ai)  = 

9  9 

and  P2<(i2<o,\<P\  (see  Definition  1). 

Fact  5. 

|^(a,)|  >  1^(6,)  |. 


(2.12) 


(2.13) 


(2.14) 


Proof.  We  prove  the  fact  for  i=l;  the  proof  for  j=2  is  virtually  identical.  From  the 
definition  of  p,  we  have  p(ai)=fl(ai)  and  p(^i)=ff(fci);  hence  —log  p(ai)  =  —log  g(ai)  and 
—log  p{bi)  =  —log  <7(61).  Therefore, 

*1 


]■[_£_(!) - Mfli)]  dx  =  0- 


a,  9 


(2.15) 


By  the  definition  of  g  the  integrand  on  the  left  side  of  (2.4.15)  is  positive  for  ai<x<Pi. 
Therefore,  the  integrand  must  be  negative  for  some  x  between  Pi  and  by  But 

min  — 2 — (x)  =  — ^ — (61) 
so  the  integrand  must  be  negative  at  by  completing  the  proof. 


Note  that 


k,  ai<x<bi 


b2<x<a2 


otherwise 


(2.16) 


which  is  non-decreasing  on  5i  =  (62,61),  non-increasing  on  S2  =  (—00,62]  L)  [61,00)  and  is  con¬ 
stant  at  its  maximum  magnitude  on  y4  =  (62,02)  U  [oi,6i). 

Fact  6  Given  two  functions  and  p^  defined  as  in  (2.15).  If 
(1)  k^  <  ku  then 

(^)  l<^fj  <  UfJ /o»- «=l-2. 

(S)  —log  Pa{x)  <  —log  Pu{x)  for  allx. 

(4)  |6,„|  >  |6,J  for  1=1,2. 

Proof.  By  Definition  1  and  Eqs.  (2.12)  and  (2.13) 

(1)  =>  (2), 

(1)  and  (2)  =>  (3), 

and  (3)  =>  (4),  completing  the  proof. 


Assume,  without  loss  of  generality,  that 

>  ^{P2) 

9  9 

Let  02  be  chosen  as  the  smallest  number  in  (P2,c)  such  that 

9  9 

Define  the  function 


^(ai)  =  /(p-ffl  0<ai<P, 


(2.17) 


(2.18) 


where  A  =  (62.02)  U  [ai,6i). 


Lemma.  If  g  and  —g'fg  are  continuous,  then  R{af)  is  a  monotonically  non- increasing  con¬ 
tinuous  function  of  a  1  with  range  (0,oo). 

Proof  See  [12] 

The  main  result  of  this  section  is  given  by  the  following  theorem. 

Theorem.  The  minimum  Fisher  Information  mixture  distribution  having  the  non-strongly 
unimodal  nominal  density  g  and  contaminant  supported  on  (62,65)  has  the  density 

/o  =  (l-f)P  +  (2.19) 

The  density  of  of  the  contaminant  is  given  by 


b,<x<a2 


(2.20) 


where 


6  +  (1-6)  j/(?  dx  +  /(?  dx 


otherwise 


(2.21) 


- 


^(oi)  =  -^(02)  . 
9  9 


and  the  constants  6,-,  j=l,2  are  found  from 

log  P(a.)  +  — (a,)  k  -  a,- 1  -  log  5(6,)  =  0  ;  (2.22) 

9 

62<i’2<02<0i<^i<6i.  If  9  and  —g'jg  are  continuous  then  a  minimum  Fisher’s  Informa¬ 
tion  distribution  in  the  form  of  Eqs.  (2.18)-(2.20)  exists  for  all  e  between  0  and  1. 

Proof.  See  [12] 


3  -  Numerical  Results 


In  §2  we  found  the  robust  detector  and  estimator  for  the  class  of  mixture  noise  distribu¬ 
tions  having  a  non-strongly  unimodal  nominal  density  g  .  Because  the  function  -g'  /g  was 
not  monotonic  we  found  it  necessary  to  limit  the  support  of  the  contaminant  to  an  interval 
about  the  mode  of  the  nominal  density.  For  distributions  satisfying  this  constraint  the  worst 
distribution  in  terms  of  Fisher’s  Information  is  given  by  Eqs.  (2.19)-(2.2l).  The  question  ^ 

which  immediately  comes  to  mind  is:  how  does  the  robust  detector  perform  if  we  violate  the 
constraint  on  the  support  of  the  contaminant?  | 

For  the  duration  of  this  section  we  will  confine  our  investigation  to  the  robust  detector; 
the  analagous  conclusions  which  can  be  drawn  for  the  estimation  problem  should  be 
apparent.  In  this  section  we  will  introduce  some  numerical  results  which  will  lend  some 
insight  into  the  question  posed  at  the  end  of  the  preceding  paragraph.  We  will  look  at  three 
different  but  related  situations.  In  the  first,  we  will  investigate  the  performance  of  several 
different  detectors  when  the  mixture  noise  conforms  to  the  restriction  on  the  support  of  the 
contaminant.  Secondly,  we  relax  this  restriction  so  that  the  interval  on  which  the  contam¬ 
inant  is  supported  is  a  superset  of  the  interval  which  was  previously  allowed.  We  then  look 
at  the  performance  of  the  same  detectors  as  in  the  first  case.  Lastly,  we  investigate  the  per¬ 
formance  of  the  robust  detector  for  several  different  distributions  which  have  support  over 
the  entire  real  line. 

All  of  the  figures  discussed  in  this  section  plot  efficacy  against  the  natural  logarithm  of 
the  mixing  constant  e.  For  simplicity,  we  assume  symmetric  noise  densities.  The  values  for 
efficacy  were  computed  using  numerical  integration. 

Figs.  1-4  plot  the  efficacies  for  the  optimal,  robust,  linear  and  sign  detectors.  The  nomi¬ 
nal  density  of  the  noise  for  these  calculated  values  is  given  by 

5  (I )  =  ^All+(|.)=1-'/"-  exp{-l/216  sinh-'(|.)f }  (3.1) 

which  is  the  expression  for  a  Johnson  density  with  parameter  6  (  see  the  discussion  of  John¬ 
son  densities  at  the  beginning  of  §2).  For  this  work  we  used  6=1  and  a=l.  The  contam¬ 
inant  density  is  given  by 


-b  <x  <b 

(3.2) 

otherwise 

where,  due  to  the  symmetry  we  have  imposed,  6  2  and  b  ^  from  (2.20)  have  become  -b  and  6  , 
respectively.  For  Figs.  1  and  3  <t^=1  and  for  Figs.  2  and  4  (T(;=10. 

As  we  know  from  §2,  the  efficacy  of  the  robust  detector  cannot  be  any  lower  for  any 
member  of  the  class  of  mixture  distributions  satisfying  the  constraint  on  the  support  of  the 
contaminant  than  it  is  for  the  worst  distribution  given  by  Eqs.  (2.19)-(2.21).  However,  this 
tells  us  nothing  about  how  well  the  detector  can  perform.  Figs.  1  and  2  show  that  for  the 
mixture  density  given  by  (3.1)  and  (3.2)  the  performance  is  not  appreciably  worse  for  the 
robust  detector  than  it  is  for  the  detector  which  is  optimal  for  those  densities.  Note  that  in 
Fig.  2  the  performance  is  slightly  worse  than  it  is  in  Fig.  1.  This  appears  to  be  due  to  the 
fact  that  for  higher  cg  more  of  the  contaminant  mass  falls  within  the  intervals  (-6  ,-o  ]  and 
jo  ,b  ).  Within  (-6  ,6 ),  the  nonlinearity  for  the  robust  detector  has  minimum  slope  and  max¬ 
imum  magnitude  on  these  two  intervab.  From  the  discussion  in  the  previous  section,  it  is 
clear  that  this  will  tend  to  reduce  the  efficacy  of  the  detector  when  compared  with  the  case 
wherein  the  contaminant  mass  is  more  concentrated  around  the  mode  of  the  nominal  density. 
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Figs.  3  and  4  plot  the  efficacy  of  the  robust  detector  for  the  noise  given  by  Eqs.  (3.1) 
and  (3.2)  along  with  that  for  the  sign  and  linear  detectors.  From  these  graphs  we  can  readily 
see  that  the  performance  of  the  robust  detector  is  far  superior  to  that  of  the  others,  espe¬ 
cially  that  of  the  linear  detector.  Both  the  linear  detector  and  the  sign  detector  have  been 
widely  used  and,  at  least  in  this  case,  it  is  apparent  that  better  results  are  readily  obtainable. 

Allowing  the  support  of  the  contaminant  to  exceed  the  the  interval  (-6  ,6  )  by  a  small 
amount  appears  to  have  little  effect  on  performance.  In  Figs.  5-8  the  contaminant  is  sup¬ 
ported  on  (-6  -  1/10,6  +  l/lO);  that  is,  it  exceeds  |  6  |  by  ten  percent  of  the  standard 
deviation  of  the  nominal,  which  in  this  case  is  0.1.  There  is  a  slight  but  largely  insignificant 
change  in  the  performance  from  the  situation  where  the  contaminant  is  confined  to  (-6  ,6  ). 

Figs.  9  and  10  compare  efficacy  for  the  robust  detector  in  two  different  mixture  noise 
environments.  The  first  noise  has  the  worst  distribution  as  described  in  §2;  the  second  has 
the  same  non-strongly  unimodal  nominal  as  the  worst  density  and  a  contaminant  which  is 
supported  on  the  entire  real  line.  If,  in  fact,  there  exists  a  distribution  which  has  minimum 
Fisher’s  Information  for  all  mixtures  with  a  given  nominal  then  it  follows  that  the  detector 
which  is  optimal  for  that  distribution  must  have  its  minimum  efficacy  for  that  same  distribu¬ 
tion;  that  is,  it  must  satisfy  the  saddlepoint  condition.  This  is  the  crux  of  the  robust  prob¬ 
lem;  to  find  a  lower  bound  on  the  performance  of  a  device  for  a  set  of  possible  circumstances 
under  which  the  device  may  be  operated.  While  we  have  not  shown  that  the  distribution 
given  by  the  theorem  of  §2  has  minimum  Fisher’s  Information  for  the  set  of  all  mixture  dis¬ 
tributions  having  a  given  non-strongly  unimodal  nominal  (we  have  done  so  for  those  with 
contaminants  having  bounded  support  on  (6o,6i))  we  can  still  find  a  lower  bound  for  the 
performance  for  distributions  having  specified  contaminants. 

In  Figs.  9  and  10  we  show  the  efficacy  for  the  robust  detector  for  the  worst  noise  distri¬ 
bution  and  for  a  noise  distribution  with  contaminant  supported  along  the  entire  real  line.  In 
both  cases  the  nominal  distribution  is  Johnson  with  6=1  and  <r=l.  The  contaminant  is  Nor¬ 
mal  in  Fig.  9  and  is  Johnson(6=l)  in  Fig.  10.  In  both  cases  the  variance  of  the  contaminant 
is  unity.  We  have  also  investigated  a  number  of  other  contaminants  including  Cauchy, 
Laplace,  Logistic,  and  Johnson  with  various  values  of  6.  The  results  are  very  similar  to  Figs. 
9  and  10.  From  the  figures  it  is  apparent  that,  for  the  set  of  contaminants  tested,  the  worst 
distribution  provides  a  lower  bound  on  performance,  as  we  would  expect. 
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